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ABSTRACT 



We prove that the classical, non-periodic Toda lattice is super-integrable. In other words, 
! we show that it possesses 2N — 1 independent constants of motion, where N is the number of 
degrees of freedom. The main ingredient of the proof is the use of some special action-angle 
(3 ■ coordinates introduced by Moser to solve the equations of motion. 
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1 Introduction 

The Toda lattice is arguably the most fundamental and basic of all finite dimensional Hamil- 
rS I tonian integrable systems. It has various intriguing connections with other parts of mathe- 
c3 I matics and physics. 

The Hamiltonian of the Toda lattice is given by 

N 1 7V-1 

H{q„ ...,qN,Pi,...,PN)=J2 7^P'+J2 ■ (1) 

1=1 ^ i=i 

This type of Hamiltonian was considered first by Morikazu Toda [1]. Equation (1) is known 
as the classical, finite, non-periodic Toda lattice to distinguish the system from the many 
and various other versions, e.g., the relativistic, quantum, infinite, periodic etc. The inte- 
grability of the system was established in 1974 independently by Flaschka [2], Henon [3] and 
Manakov [4] . The original Toda lattice can be viewed as a discrete version of the Korteweg- 
de Vries equation. It is called a lattice as in atomic lattice since interatomic interaction was 
studied. This system also appears in Cosmology. It appears also in the work of Seiberg and 
Witten on supersymmetric Yang-Mills theories and it has applications in analog comput- 
ing and numerical computation of eigenvalues. But the Toda lattice is mainly a theoretical 
mathematical model which is important due to the rich mathematical structure encoded in 
it. 

Hamilton's equations become 



The system is integrable. One can find a set of independent functions {Hi, . . . , H^} whicli 
are constants of motion for Hamilton's equations. To determine the constants of motion, 
one uses Flaschka's transformation: 



(2) 



Then 
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These equations can be written as a Lax pair L — [S, L], where L is the Jacobi matrix 
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This is an example of an isospectral deformation; the entries of L vary over time but the 
eigenvalues remain constant. It follows that the functions Hj — jtrL-' are constants of 
motion. The mappings (2) is a projection from R^^ to R^^~^ and it is clearly not one-to- 
one. Replacing the vector q by q + c gives the same image. 

The supper-integrability of this type of systems should be expected due to their dispersive 
asymptotic behavior. However, the construction of integrals is not typically a trivial task. 
In the case of the open Toda lattice, asymptotically the particles become free as time goes 
to infinity with asymptotic momenta being the eigenvalues of the Lax matrix. Therefore, 
the system behaves asympotically like a system of free particles which is sup er-integr able. 
The super-integrability of the Toda lattice ioY N — 2 was estabhshed in [5] . The additional 
integral in [5] was obtained using Noether's theorem. We give the formula for the additional 
integral. 



In the case of two degrees of freedom the potential is simply 

V{qi,q2) = e'^'-'^- , 
and the procedure of Noether produces the following three integrals: 

h = 7= exp ^ Ji ■ . (4) 

P1-P2- vJi V P1 + P2J 

Note that H = Hf + | Ji and that the function G = -^^^ which appears in the exponent of 
Ji is a time function, i.e., it satisfies {G,H} = 1. We note that the integral Ii remains an 
integral if we add a real constant to qi + q2, an observation that will be important later on. 
The existence of the integral Ii shows that the two degrees of freedom Toda lattice is super- 
integrable with three integrals of motion {Hi, Ji,Ii}. As we will see, the complicated integral 
Ii has a simple expression if one uses Moser's coordinates. 



2 Moser's solution of the Toda lattice 



Moser's beautiful solution of the open Toda lattice uses the Weyl function /(A) and an old 
(19th century) method of Stieltjes which connects the continued fraction of /(A) with its 
partial fraction expansion. The key ingredient is the map which takes the (a, b) phase space 
of tridiagonal Jacobi matrices to a new space of variables (Aj, rj) where Aj is an eigenvalue 
of the Jacobi matrix and r, is related to the residue of some rational functions that appear 
in the solution of the equations. We present a brief outline of Moser's construction. 
Moser in [6] introduced the resolvent 

R{\) = {XI -ly^ , 

and defined the Weyl function 

/(A) = (i?(A)eiv, ejv) , 

where bat = (0, 0, . . . , 0, 1). 

The function /(A) has a simple pole at A = Aj. For the purpose of this paper, and this is 
also observed by Moser, we define / by the formula 

/ N 2 \ N 

The differential equations in the variables (A, r) take a particularly simple form: 

A, = 

h = -Aj n . (6) 
These equations show that (Aj, logr^) are action-angle variables for the Toda lattice. 



(5) 



The variables of, hi may be expressed as rational functions of Aj and rj using a continued 
fraction expansion of /(A) which dates back to Sticltjcs. Since the computation of the 
continued fraction from the partial fraction expansion is a rational process the solution is 
expressed as a rational function of the variables (Aj, rj). The procedure is as follows: 
The Rnn element of the resolvent, as defined previously, takes the following continued 
fraction representation: 



/(A) = 
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The function / (A) has N simple poles at the eigenvalues of the Lax pair matrix L. Therefore, 
its partial fraction expansion has the form: 



/(A) = E 



N ^2 



N 



(8) 
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where the residue of /(A) at A = Aj is rf / {J2i^i i^f) ■ Stieltjes described a procedure that 
allows one to express and 6j in terms of Ai, . . . , A at and ri, . . . , r^r. We briefly describe 
the method. We expand the partial fraction expansion of /(A) as given in (8) in a series of 
powers of j. We obtain. 
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The coefficient of A^"*"^ is denoted by cj and equals. 
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The formulas of Stieltjes involve certain i x i determinants which we now define. 
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The formulas that give the relation between the variables (a, 6) and (r. A) are. 
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^N-i 
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= 1 AT 



where = 1, -Bo = 1, -B-i = 0. 
For example, in the case N — 2 



^1 — Co, A2 — CqC2 — C^, Bi — Ci, B2 — C1C3 — C2 



and therefore 



Thus, 



A2 -Bo 



a 



r^i(A2-Ai)^ 

^ !i^l±4^ (9) 
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One can check that the differential equations fi — —riXi, for i = 1, 2 correspond via trans- 
formation (9) to the A2 Toda equations 



di = 01(62 - bi) 
hi — 2af 
62 = -2ai . 

As Moser notes, it is not too hard to obtain exphcit expressions for = 3 but the general 
case is quite comphcated. With the exception of A" = 2 the a, are not rational functions of 
(Ai,ri) but the are. In general one can express (at least in theory) the functions Oj, bi in 
terms of rj, Aj. Again, the function is not one-to-one. In fact, if one replaces the vector r 
with ar the result is the same. 

Finally, we comment on the Poisson brackets in the new coordinates. The multi-hamiltonian 
structure of the Toda lattice was developed in [7] and [8] using master symmetries. The 

analogous results in (A, r) coordinates are due to Fcybusovich and Gekhtman [9]. The Poisson 
brackets project onto some rational brackets in the space of Weyl functions and in particular, 
the Lie-Poisson bracket of the Toda lattice corresponds to the Atiyah-Hitchin bracket [10]. 
In general, one constructs a sequence of Poisson brackets on the space (Aj,rj) whose image 
under the inverse spectral transform corresponds to the standard Toda hierarchy. A rational 
function of the form is determined uniquely by the distinct eigenvalues of p(A), Ai, . . . , A„ 

and values of q at these roots. The residue is equal to and therefore we may choose 

Ai,...,A„,g(Ai),...,g(A„) 

as global coordinates on the space of rational functions (of the form | with p having simple 
roots and q,p coprime). We have to remark that the image of the Moser map is a much 
larger set. 

The kth Poisson bracket is defined by 



{A„ g(A,)} = -Xhi\) 
{q{^^), 5(A,)} = {A„ A,} = 



The initial Poisson bracket under the Moser map is given exphcitly by 



{A,, A,} = 
{ri,rj} = 

{Xi, Tj} = SijTj i,j ^1,...,N . (10) 

Similarly, the quadratic Toda bracket, corresponds to a bracket with only non-zero terms 
{Kn} = XiVi. 

1 ^ 

The Hamiltonian function in the new coordinates is H2 = - ^X^. In other words, taking 

2 i=i 

H2 as the Hamiltonian and using bracket (10) gives equations (6). 

3 The Toda lattice is super— integrable 

We now come to the main result of this paper. We define 

= e^-+\ i = l,...,7V-l, (11) 

where 

It is easily shown, using equation (6) that — - = 0, for j = 1, . . . , N — 1 and thus the 

at 

functions Ij are constants of motion. 

The functions Hi — X\ + Xl + . . . + X]^ and Ij, i — 1, . . . ,N , j — 1, . . . ,N — 1 are functionally 
independent. In fact, the Jacobian (2A^ — 1) x 2N matrix of the functions Hi and Ij has a 
(2A^ — 1) X (2A^ — 1) subdeterminant, dN+i-, which is obtained by deleting the (A'"+ l)-column 
and is not identically zero. A simple calculation gives 

d^+i = -2^-^7V ^^^e^^.- n (A. -A,). 

'jv-2'jv-i'Ar 1 l<i<j<N 

Since the eigenvalues of real Jacobi matrices are distinct, the functions Hi and Ij are inde- 
pendent. We summarize the results in the following: 

Theorem 1 The Toda lattice with N degrees of freedom posesses 2N — 1 independent con- 
stants of motion, Hi, i — 1, . . . , N , Ij, j — 1, . . . , N — 1, and is therefore super-integrable. 

Remark 1 

It is clear that the functions n — 1, . . . ,N are in involution. Moreover it can be shown 
that {/j, Ij} — 0, i,j — 1,...,N — 1. In addition, lor n — 1, . . . ,N, j — 1,...,N — 1 



where c„ = 1 for n = 2, . . . , A'", ci = and 



Ej = E - (E ^0 - 2) - 2A,+iA, w{n - 3). 



The sums are taken over all i from 1 to where i ^ j, j + 1 The function w{n) symbolizes 
the full homogeneous polynomial in Xj and A^+i that have total weight equal to n. For 
instance, w{n) = 0, n e — Z+, w{Q) = 1, w{l) = Xj + A^+i, w{2) = A| + A^_,_i + AjAj+i, etc. 
One can, of course, use the quadratic Toda bracket in {Xi,ri) coordinates. We must then 
take TrL = Ai + . . . + A^v as Hamiltonian. However, in this bracket the Hi, Ij do not form a 
finite dimensional algebra. 

Remark 2 

We clearly have {H2, Ij} = 0, j — 1, . . . , N—1, since H2 is the Hamiltonian and the functions 

Ij are constants of motion. 

We define the sets 5*1 = {Hi, . . . , Hjsi} and ^2 = {H2, h, ■ ■ ■ , In-i}- Then if /, g( e =^ 
{f:9} — ^-iid if £ 'S'2 =^ {f:9} — 0- III other words the sets Si and S2 are both 
maximal sets of integrals in involution. We therefore have two different sets demonstrating 
the complete integrability of the Toda lattice. 

Remark 3 

We finally would like to comment on how the integrals Ij were guessed: The complicated 
integral (4) at the end of the introduction is quite simple in Moser's coordinates. For example, 
^/Jl is simply equal to 2(A2 — Ai) and the expression 

P1-P2 + yji 

P1-P2- ^\ 

reduces to — . The exponent is simphfied as follows. On the one hand, 

(?i + 52)' = Pi + P2 = -2(61 + 62) = -2(Ai + A2) . 

On the other hand from rj = — nAj one obtains that (Inrj) = —A,. Therefore, the function 
gi + q2 and 21nrir2 differ only by a real constant. In other words, up to a constant, the 
exponent is simply 

2(Ai - A2) , , . 

^^TTaT'"^'^'^^' 

and up to a sign difference, the integral (4) is precisely the same as the one in (11). 
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